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, Abstract 

C**) ' Using ODE techniques we prove the existence of large classes of initial 

[ data satisfying the constraints for the spherically symmetric Einstein- 

■ Vlasov-Maxwell system. These include data for which the ratio of total 

charge to total mass is arbitrarily large. 
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1 Introduction 



The global dynamical behavior of self-gravitating matter is a subject of cen- 
tral importance in general relativity. A form of matter which has particularly 
■ nice mathematical properties is collisionless matter, described by the Vlasov 

equation. It has the advantage that it lacks the tendency observed in certain 
other models, such as perfect fluids, that solutions of the equations of motion 
of the matter lose differentiability after a finite time. These singularities of the 
mathematical model form an obstacle to further analysis and prevent the study 
of the global dynamical properties of the solutions. Collisionless matter is free 
from these difficulties and there is a growing literature on global properties of 
solutions of the Einstein- Vlasov system PP , [S| ■ 

In [7], the authors prove the global existence of asymptotically flat solu- 
tions of the spherically symmetric Einstein- Vlasov system, with small initial 
data. That study concerns uncharged particles. We consider, under the same 
assumption of spherical symmetry, the case where the particles are charged. To 
describe the full physical situation, we must then couple the previous system 
to the Maxwell equations that determine the electromagnetic field created by 
the fast moving charged particles. As will be seen below, that reduces, in the 
spherically symmetric case, to its electric part. 
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It is appropriate at this point to examine the motivation for considering this 
particular problem. We are not aware that it has any direct astrophysical appli- 
cations. There are, however, two reasons why the problem is interesting. The 
first is that it extends the knowledge of the Cauchy problem for systems involv- 
ing the Vlasov equation and it will be seen that it gives rise to new mathematical 
features compared to those cases studied up to now. The second is connected 
with the fact that it would be desirable to extend the work of [7] beyond spher- 
ical symmetry. In particular, it would be desirable from a physical point of 
view to include the phenomenon of rotation. Unfortunately, presently available 
techniques do not suffice to get away from spherical symmetry. In this situa- 
tion it is possible to attempt to obtain further intuition by using the analogy 
between angular momentum and charge, summed up in John Wheeler's state- 
ment, 'Charge is a poor man's angular momentum'. Thus we study spherical 
systems with charge in the hope that this will give us insight into non-spherical 
systems without charge. This strategy has recently been pursued in the case of 
a scalar field as matter model, with interesting results [2], 0. 

Concerning the Cauchy problem in general relativity, it is well known that 
in addition to the Einstein evolution equations there are constraint equations 
which have to be solved. (See e.g. @|.) It is only of interest to consider the 
problem of evolution, if the problem of constraints on the initial data can be 
solved. In our specific case, we are led to a difficulty in solving the constraints 
on the initial data, that has not previously been considered in the literature. 
Let us first recall the situation in [7j before seeing how it changes in the case 
of charged particles. In 0, using the assumption of spherical symmetry, the 
authors look for two metric functions A and /i, that depend only on the time 
coordinate t and on the radial coordinate r, and for the distribution function / 
of the uncharged particles that depends on t, r and on the 3-velocity v of the 
particles; the metric functions A, fj, are subject to the Einstein equations with 
sources generated by the distribution function / of the collisionless uncharged 
particles which is itself subject to the Vlasov equation. They show that the 
Einstein equations to determine the unknown metric functions A and fi, turn 
out to be two first order ODE in the radial variable r, coupled to the Vlasov 
equation in /. Putting t = in the Einstein equations yields two constraint 
equations that link the three initial data that are the two initial data for A(t, r) 

and n(t, r) denoted A(r) and fi(r) respectively, and the initial datum for f(t, r, v) 

° o 

denoted by f(r, v). The equations to be solved are two first order ODE for A, /i. 
The mass function m(t, r) is defined as an integral of / over the hypersurfaces of 
constant t. The exact definition is given in the next section. Provided the initial 

mass function m(r) = m(0,r) is everwhere less than r/2 the function A can be 
determined from the relation 1 — 2m/ V = e~ 2A . Once this has been done ^ can 

° o 

be determined by integration. In this way functions A and /i can be determined 

o 

in a straightforward way when / is given, provided the inequality 2m /r < 1 is 
satisfied on the initial hypersurface. Thus a simple parametrization of the set 
of initial data satisfying the constraints is obtained. They can be constructed 
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from a non-negative function / which is required to satisfy one inequality. 

In the case of charged particles, the problem of constraints on the initial data 
is not so easy. We consider the case of a spherically symmetric electric field E of 
the form E(t, r) = e(t, r)£, where e(t, r) is an unknown scalar function and r the 

position vector in R 3 . We denote by e the initial data for e(t,r). The Maxwell 
equations imply a constraint equation on the initial data, that is a first order 
ODE in the radial variable r. We now have to face the problem of constraints 

°o°o • r \ °i o 

between the four initial data A, fi, /, e. If the equations for A, /, and e can 
be solved then fj, can be determined by a simple integration as in the case of 
uncharged particles. On the other hand the problem of determining the first 
three quantities is more difficult in the charged case due to the coupling of the 
gravitational and electromagnetic constraints. Any solution of the constraints 
must satisfy the condition 2m/ r < 1 on the initial hypersurface. However in 
the charged case m includes a contribution due to the energy density of the 
electric field and this depends on e and A. As a consequence, in contrast to the 
uncharged case, this condition cannot be expressed as an inequality on / alone. 
Moreover, the equation for e also contains A and so it is not possible to solve 
the electromagnetic constraint separately from the gravitational one. There is 

o 

no alternative but to solve the two constraints together. When / is prescribed 
the constraints define a system of two coupled nonlinear ordinary differential 

equations for A and e. These equations are singular at r = 0. 

Solving the ordinary differential equations arising from the constraints is a 
task which must necessarily be mastered before undertaking the study of the 
problem of evolution for the Einstein- Vlasov-Maxwell system. Since the equa- 
tions are singular they cannot be handled by standard results of ODE theory 
alone. In this paper, we prove, using statements developed in |U] on singular 
ODE, that, under smallness assumptions on the prescribed initial data function 

o 

/, a corresponding global solution of the constraints exists. This provides a way 

of showing the existence of initial data A, //, e, in the case of an asymptotically 
flat spacetime with a regular center. In other words, we prove that the results 
obtained for the constraint equations in jj] that focuses on the case of uncharged 
particles, extend to the case of charged particles, in the sense that the initial 

o 

datum / can still be considered as the only arbitrary initial datum. 

Given the complicated nature of the constraint equations it seems impractible 

o 

to describe the most general class of functions / for which the constraints can 
be solved. What we have been able to do is to describe two large classes of 

o 

functions / for which it is possible. The first case treated is a rather obvious 
one to try. The charge q of a particle is a free parameter in the equations. In 
the case q = the equations reduce to the uncharged situation analysed in [Jj. 

o 

A perturbation argument then allows the case with fixed / and small q to be 
treated. A more exotic class of initial data is obtained by a rigorous perturbation 
argument in a regime where the density of particles tends to zero and the charge 
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per particle tends to infinity. This allows the construction of initial data which 
coincide outside a compact set with data for the Reissner-Nordstrom solution 
with arbitrarily large charge to mass ratio. The restriction of these data to the 
region outside a large radius does not allow an electrovacuum interior which has 
a regular centre or corresponds to a black hole but, as this construction shows, 
does allow a regular (dynamical) interior with charged matter. 

The paper is organized as follows. In Sect. 2, we recall the spherically sym- 
metric Einstein- Vlasov-Maxwell system, from which we deduce the constraint 
equations. In Sect. 3, we recall the main result of [5j which we use and we prove 
the existence of large classes of solutions of the constraint equations for the 
Einstein- Vlasov-Maxwell system. The derivation of some facts required in Sect. 
3 concerning regions where there are no particles is contained in an appendix. 



2 Formulation of the problem 

We consider fast moving collisionless particles with unit mass and charge q. The 
basic spacetime is (M 4 , <?), with g a Lorcntzian metric with signature (— , +, +, +). 
In what follows, we assume that Greek indices run from to 3 and Latin 
indices from 1 to 3, unless otherwise specified. We also adopt the Einstein 
summation convention. The metric g reads locally, in cartesian coordinates 
{x a ) = (x°,x l ) = {t,x): 

ds 2 = g aP dx a <g> dx p (2.1) 

The assumption of spherical symmetry means that we can take g of the following 
form (Schwarzschild coordinates) 

ds 2 = -e 2tl dt 2 + e 2X dr 2 + r 2 (d9 2 + sin 2 6dip 2 ) (2.2) 

where fi = fi(t,r); A = A(t,r); t G K; r G [0,+oo[; 6 G [0,tt]; <p G [0,2tt]. The 
Einstein- Vlasov-Maxwell system reads: 

RaP ~ \9 a pR - 87T(T a/3 (/) + T af){F)) (2.3) 

£x(F)f - (2.4) 
V a F af) = J?; V a F^ + V fi F ia + V 7 F Q/3 = (2.5) 



with: 



T a0 (f) = - [ Pa ppfu p ; T a p{F) = - 9 JfF lv F^ + F 01 F a ~< 
Jr 3 4 

J P {f){x) = q [ p P f{x,p)uj p , uj p =| g \ * dp dp dp , p = g oP° , 
Jm 3 Po 

X a (F) - (p a ,-r^ P V - qp^Fp ), 
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where denote the Christoffel symbols. Here, x = (x a ) is the position andp = 
(p a ) is the 4-momentum of particles. In the expressions above, / stands for the 
distribution function of the charged particles, F stands for the electromagnetic 
field created by the charged particles. Here l|2.3[l are the Einstein equations for 
the metric tensor g = {g a p) with sources generated by both / and F, that appear 
in the stress-energy tensor T a p + T a p. Equation i|2.4[l is the Vlasov equation for 
the distribution function / of the collisionless particles and (|2.5(l are the Maxwell 
equations for the electromagnetic field F, with source (current) generated by / 
through J = J(f). One verifies that the conservation laws V Q ,(T Q ' 3 + r a ^) = 
hold if / satisfies the Vlasov equation. 

By the assumption of spherical symmetry, we can take g in the form 1|2.2|) . 
One shows, using the Maxwell equation that F reduces to its electric part. We 
take it in the form E = (E a ) with E° — 0, E l — e(i,r)^r, and then, a straight- 
forward calculation shows that: 



too = -e 2 ^e 2 (t,r); r 0l = 

where 5y is the Kronecker symbol. 

These relations and results of 7 show that the spherically symmetric Einstein- 
Vlasov-Maxwell system implies the following first order ODE system in A, fi, /, 
e: 

e" 2A (2rA' - 1) + 1 = %irr 2 p (2.6) 
e - 2X (2rfi' + 1) - 1 = 8vrr 2 p (2.7) 



at y/l + v 2 ox r r ov 

(2.8) 

— (r 2 e x e{t,r))=qr 2 e x f(t,x,v)dv (2.9) 



where A' = S*; A = f£ and: 



dr 
at 



p(t,x) = [ f(t, x, v)V~l + v 2 dv + l e 2X(t ^e 2 (t, x) 

JR3 2 



(2.10) 



p{t,x)= f (^) 2 /(^i,,)^L=-ie 2 ^) e 2 (^) (2.11) 

Here <|2.6[1 - <|2.7[) are the Einstein equations for A and p, l|2.8l) is the Vlasov 
equation for / and l|2.9|l is the Maxwell equation for e. Here i and v belong 
to M 3 , r := \x\, x ■ v denotes the usual scalar product of vectors in M 3 , and 
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v 2 := v ■ v. The distribution function / is assumed to be invariant under 
simultaneous rotations of x and v, hence p and p can be regarded as functions 
of t and r. It is assumed that f(t) has compact support for each fixed t. We 
are interested in asymptotically flat spacetimes which leads to imposing the 
boundary condition that: 

lim fx{t,r) = (2.12) 

r — >oo 

They should also have a regular centre which means that in addition to A, \x and 
e being smooth functions oft and r, including at r = 0, the boundary condition 
X(t, 0) = should be satisfied for all t. Note that the condition that A(t, r) 
tends to zero as r — > oo, which is part of asymptotic flatness, follows from the 
field equations and the fact that /(£) has compact support. This is because in 
the region where / vanishes the general solution of 12. 6|) satisfies the condition 
A(r) = 0{r) as r — > oo. 

Let m(r) = 4-7T J Q r s 2 p(s)ds. This is the mass function referred to in the 
introduction. Its limit M as r — > oo is the total or ADM mass of the system. 
The function n = J R3 fdv is the number density of particles and nq the charge 
density. The total charge of the system is given by Q — Anq J °° s 2 e x ^n(s)ds. 

Now, define the initial data by: 

/(0, x, v) = f(x, v); A(0, x) = X(x) = A(r); 1(0, x) - e(x) - e(r) (2.13) 

o 

with / S Cg being a continuously differentiable function with compact support, 
which is nonnegative and spherically symmetric, i.e 

MA e 50(3), V(£,«) e R 6 , f(Ax,Av) = f(x,v). 

We obtain the constraint equations on the initial data by taking l|2.(j|) , l|2.7|l and 
fHfy for t = 0, that give: 

e- 2X (2r°X - 1) + 1 = 87rr 2 \ y/l + v 2 f(r, v)dv + Avr 2 e 2A t (2.14) 
±(r 2 e l °e) = qr 2 e l f °f(r, v)dv = J(f) (2.15) 

* JR3 

e- 2l (2r;i' + 1) - 1 = 8nr 2 [ (— ) 2 /(r, ») - = - 4^r 2 e 2 ° x t (2.16) 

Jr3 r VI + v 2 

o o 

We observe that, if / is given and if we can solve 1)2. 14(1 - 1)2.15(1 for A and e 

o ' ... 

then H2.1(j|) determines at once /i . Using the boundary condition ()2.12|) then 

determines p. So we can concentrate on H2.14|) - (|2.15(l . In what follows, we fix 

o 

/ in H2.14|) - ()2.15|) and we look for a unique global asymptotically flat solution 

o 

(A, e) of the system ^2.14^ - ^2.15^ above with regular centre. Note that, using 
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the compact support of / and the equations l|2.14|) and H2.15|l . it follows that 

A and e tend to zero as r — > oo. It also follows from H2.14(l and 1)2. 15[) and the 

°' 

regularity of the solution that A (0) = and e(0) = 0. 

3 Existence of global solutions of the constraints 

In this section the existence of global solutions of the equations (|2.14|) and l|2.15|l 
will be proved. Let us state first of all the following result of |5] on which our 
global existence theorem relies: 

Theorem 3.1 Let V be a finite- dimensional real vector space, N : V — > V a 
linear mapping, G : / x V — > V a smooth (i.e C°° ) mapping and g : I —> V a 
smooth mapping, where I is an open interval in R containing zero. Consider 
the equation 

s^ + Nf = sG(s,f(s)) + g(s) (3.1) 

for a function f defined on a neighborhood of in I and taking values in V. 
Suppose that each eigenvalue of N has a positive real part. Then there exists an 
open interval J with € J C / and a unique bounded C 1 function f on J \ {0} 
satisfying Moreover f extends to a C°° solution of \ H. 1\) on J . If N , G 

and g depend smoothly on a parameter z and if the eigenvalues of N are distinct 
then the solution also depends smoothly on z. 

Proof. See Theorem 1 in 0, p.989. 

Remark 1 The assumption that N has distinct eigenvalues is to ensure that 
N can be reduced to diagonal form by a similarity transformation depending 
smoothly on z. In particular, Theorem 3.1 applies if N is already a diagonal 
matrix. 

o 

Theorem 3.2 (Local existence) Let f 6 C°°(R 6 ) be nonnegative, compactly 
supported and spherically symmetric. Then, the equations {2.1$ and \2.15)) 

have a unique local and regular solution (A, e) defined on some interval [0,i?], 
R > 0. The solution depends smoothly on the parameter q. 

o 

Proof: Let / 6 C°°(M. 6 ) be nonnegative, compactly supported and spherically 
symmetric. By a regular solution we mean one which is smooth and for which 

o o 

A = at r = 0. It follows that for any regular solution A can be written in the 
form: 

A(r) = rL{r) (3.2) 

o' 

for a smooth function L(r). Equation (|3.2|l implies A = L + rL' and l|2.14|l - 
(|2.15|l can be written: 

rL' + L = — (l-e 2rL )+4irre 2rL { [ y/l + v 2 °f(r, v)dv + -e 2rL t) (3.3) 
2r Jr3 2 
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re + 2e = -re(L + rL') + rq l f(r,v)dv (3.4) 

Jn 3 

The function e 2x — 1 — 2x vanishes to hrst order at the origin and hence e 2x — 
1 — 2x = x 2 Fq(x) for a smooth function Fq. Hence the equation for L can be 
rewritten in the form 

rL' + L = -L+ r -L 2 F {rL) + ^re 2rL { [ ^Jl + v 2 °fdv + \e 2rL t) 
2 Jr3 2 

Thus 

rL' + 2L = rGi(r,L,e,f) (3.5) 

where 

Gi(r, L, e, /) = )-L 2 F (rL) + 4vre 2 '' L / y/l + v 2 f(r, v)dv + 27re irL t 

and (|3.4|l reads, given (|3.5|) : 

re + 2e = — re(— L + rG\(r, L, e, /)) + rq I f(r,v)dv. 

Jr 3 

Hence: 

re + 2°e = rG 2 {r,L°eJ) (3.6) 

where 

G 2 (r,L,e, f) = Le- reGi(r, L,e, /) + g / f(r,v)dv 

Jn 3 

^ J and $ - ]/ ^ 
and (|2.15|) can be written: 



Setting G = ( G J and $ = ( o ) and using H3.5H3.6j l. tne equations l|2.14jl 



r — + 2$ = rG(r,$(r)) (3.7) 
ar 

We apply Theorem 3.1 with V = tt 2 , TV = ^ ^ , 7V$ = 2$ to and, 

since G clearly depends smoothly on q, obtain the desired result. Thus Theorem 
3.2 is proved. 

o 

Theorem 3.3 (Global existence, low charge) Let f E G°°(M 6 ) be nonneg- 
ative, compactly supported and spherically symmetric with 



8ir s 2 ( y/l + v 2 f(s,v)dv)ds < r. (3.8) 



Then, for q small enough, the equations {2.14V and $2.15\) have a unique global 
and regular solution (A, e) defined on [0, +oo[ that satisfies the boundary condi- 
tion A(0) = 0. 
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Proof Let / £ C°°(M. 6 ) be nonnegative, compactly supported and spherically 

o 

symmetric. We assume that / is fixed and satisfies l|3.8(l . By Theorem 3.2, the 
equations H2.14(l and H2.15[l have a unique local regular solution on some interval 

o 

[0, .R], R > 0. Again, Theorem 3.1 shows that, for fixed /, there exists E > 0, 
such that for q £ [-E, E], R can be chosen uniformly and the solution on [0, R] 

o 

depends continuously on the parameter q. Now, for fixed / and q, the solution 

o 

has a right maximal interval of existence [0,i2*[, i?,* = R*(f,q). We have to 
show that i?* = +oo. In fact, the second term in the r.h.s of (|2.14|l vanishes for 
q = 0, as one can see by integrating Ij2.15|l over [0, r], r > 0. It follows that for 
q = 0, l|2.14|l and (|2.15|) have a global solution under the sole assumption 1|3.8JI 

o 

on /. Then by the stability theorem for ODE, for every R > 0, there exists a 
number E > 0, such that, for every q £ [-E, E], the system (13.711 has a solution 
$ s that exists on [0,R] (see Theorem 4, p. 92 in fj). Thus i?* > R. Now, we 

o o 

can choose R large so that supp/ C [0, R] x M 3 , i.c f(r, v) = for r > R. If i?, 
is the radius of the support of the distribution function then R may be chosen 
to be bigger than m(R ) + Q 2 /(8nR ) for all q in the interval [E, E]. Hence 
by the lemma of the appendix the solution extends to one which is global and 
regular. This completes the proof of the theorem. 

Theorem 3.4 (Global existence, high charge) Let f £ C°°(IR 6 ) be non- 
negative, compactly supported and spherically symmetric. Then, for q large 

enough, the equations \B. J^| ) and ^2.15)) have a global and regular solution (A, e) 

o o 

defined on [0, +oo[ that satisfies the boundary condition A(0) = for which f is 
a constant multiple of f . Moreover the charge to mass ratio Q/M of the solution 
can be made as large as desired. 

Proof We assume that / is fixed as in the assumptions of the theorem We set: 

a = q- 1 ; f = aT*/; e = a _( * -1) e 
for some integer k > 2. Then Q2.14[l - l)2.15[l can be written as: 

e~ 2 ^(2rA - 1) + 1 = 8nr 2 (a k f f(r, v)dv + \e 2 ° x a^ k ^e 2 ) (3.9) 

2e + reX + re' = -r [ f(r,v)dv. (3.10) 

Introducing a variable L as defined in l|3.2(l puts these equations into a form 
closely analogous to that obtained in the proof of the last theorem. In fact 
the left hand side has the same form as in that case. All that is changed is 
the form of the nonlinear terms on the right hand side. The equations depend 
on a as a parameter in a way which is smooth in a neighbourhood of a = 0. 

o 

For a = the function A vanishes identically while the equation for e becomes 
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linear and has a global regular solution. From this point on we can argue just 
as in the proof of Theorem 3.3 to conclude that for a sufficiently small there is a 
unique global regular solution of these equations. Here we must use the fact that 
rri(Ro) + Q 2 /(8nR ) is bounded independently of a for a small. The assumption 
that a is small corresponds to q being large. The distribution function belonging 

o 

to the solution is obtained from / by a constant rescaling. The charge to mass 
ratio of the solution is proportional to a~ k and thus tends to infinity as a tends 
to zero. This completes the proof. 

Remark 2 The solution in the exterior region is part of the Reissner-Nordstrom 
solution JhTj. 

Remark 3 Our motivation in proving these theorems was to construct initial 
data for the Einstein- Vlasov-Maxwell system. The same arguments apply with 
other kinds of charged matter such as a charged fluid as sources for the Einstein 
equations. 

Acknowledgments: One of us (ADR) thanks Curt Cutler for a helpful sugges- 
tion. This work was supported by a research grant from the VolkswagcnStiftung, 
Federal Republic of Germany. 

A Appendix 

This appendix contains an analysis of exterior regions free of particles in initial 

o 

data sets for the Einstein- Vlasov-Maxwell system. Let / be an initial distri- 
bution function of compact support. Let Ro be the radius of its support in 

o 

space so that / vanishes for all r > Ro- In this appendix we are only concerned 
with quantities on the initial hypersurface and so we will drop the label zero 
indicating the restriction of spacetime quantities to the initial hypersurface. 

Lemma 1 Consider a solution of the constraint equations for the spherically 
symmetric Einstein-Vlasov-Maxwell system defined for § < r < R\ and having a 
regular centre. Suppose that radius Ro of the support of the distribution function 
f is less than R\. Let M = m(i?o) + Q 2 1 '(871 'Ro) . Then if Ri > 2M the given 
solution extends to a unique solution of the constraints defined for all r £ [0, 00 [ 
which is asymptotically flat and has f = for R > R . 

Proof Integrating the constraint equation Ij2.15|l gives 

r 2 e x{r) e{r)=q f s 2 e x{s) ( f(s,v)dv (A.l) 
Jo Js. 3 

For r > Rq the upper limit r in the integral can be replaced by Ro or infinity 
without changing the value of the expression. It is equal to Q/Air where Q is 
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the total charge of the system defined in Sect. 2. For r > Rq the function / 
vanishes and the mass function m defined in Sect. 2 satisfies 

2-7T 

m' = —(Q/4ir) 2 (A.2) 

It follows that M(r) = m(r) + Q 2 /(%irr) is independent of r. If the solution 
exists globally in r and is asymptotically flat then taking the limit r — > oo shows 
that M is equal to the ADM mass M. In any case M is positive and it follows 
that in the exterior region m = M — Q 2 /(8irr). In order to determine whether 
a solution can be extended to larger values of the radius it is enough to ensure 
that 1 — 2M/r + Q 2 /(Airr 2 ) remains positive. For in that case we can define A 
by means of the relation 

e" 2A = 1 - 2M/r + Q 2 /{Ai:r 2 ) (A.3) 

Note that linv^oo A(r) = 0. Once this has been done we can define fi to be 
equal to —A and e(r) = r~ 2 e~ x {Q / Air) in the external region and this gives 
the unique solution satisfying the correct boundary conditions. If r > 2M 
then 1 — 2M/r + Q 2 / (&irr 2 ) is automatically positive and the desired result is 
obtained. 
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